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ABSTRACT 

We have simnlated the Expanded Owens Valley Solar Array (EOVSA) radio 
images generated at mnltiple freqnencies from a model solar active region, em¬ 
bedded in a realistic solar disk model, and explored the resnlting datacnbe for 
different spectral analysis schemes to evalnate the potential for realizing one of 
EOVSA’s most important scientific goals—coronal magnetography. In this paper, 
we focns on modeling the gyroresonance and free-free emission from an on-disk 
solar active region model with realistic complexities in electron density, temper- 
atnre and magnetic field distribntion. We compare the magnetic field parame¬ 
ters extrapolated from the image datacnbe along each line of sight after folding 
throngh the EOVSA instrnmental profile with the original (nnfolded) parameters 
used in the model. We find that even the most easily automated, image-based 
analysis approach (Level 0) provides reasonable quantitative results, although 
they are affected by systematic effects due to hnite sampling in the Fourier (uv) 
plane. Finally, we note the potential for errors due to misidentified harmonics of 
the gyrofrequency, and discuss the prospects for applying a more sophisticated 
spectrally-based analysis scheme (Level 1) to resolve the issue in cases where 
improved uv coverage and spatial resolution are available. 

Subject headings: Sun: corona — Sun: magnetic fields — Sun: radio emission 


1. Introduction 

The magnetic held in the solar atmosphere plays a critical role in determining the 
plasma structure, storage of free magnetic energy, and its release in various forms of solar 
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activity, such as flares and coronal mass ejections (CMEs). In order to interpret such solar 
phenomena, the calculation of free magnetic energy and its storage-release process requires 
quantitative knowledge of the three-dimensional (3D) structure of the corona in or near the 
span of eruption. 


One way to derive the coronal held strength is through nonlinear force-free held (NLFFF) 
extrapolations from the measured photospheric boundary data. To further verify the model, 
morphological tests between the simulated held lines and 2D coronal brightness maps ob¬ 
served in EUV and soft X-rays are performed. However, these extrapolations are found to be 
imperfect when applied to real solar data: different algorithms demonstrate various coronal 
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which may affect the boundary become critical to the success of NLFFF modeling in real 
applications: (1) measurement uncertainties, such as noise levels and the 180-degree ambi¬ 
guity, (2) measurement limitations, specifically deriving the vector magnetic held over a held 
of view large enough to accommodate most of the field-line connectivity within the active 
region and its surroundings, and insufficient spatial resolution to recover the distribution of 
current densities on the boundary, and (3) poor knowledge of an appropriate procedure for 
transforming the non-force-free photospheric magnetic fields into an approximately force-free 
boundary. Fven if some of these issues can be addressed, NLFFF extrapolation relies heavily 
on the assumption that the coronal field is force-free, which is not always true in the case 
of fast eruptive activities. Moreover, morphological tests are constrained by the 2D obser¬ 
vational clues to the coronal structure, where the underlying field lines may not necessarily 
follow the observed loop structure in brightness. Therefore, it is essential to develop direct 
measurement techniques in conjunction with extrapolations to better determine the coronal 
magnetic field. 


Unfortunately, it is very difficult to directly measure the coronal magnetic field through 
conventional Zeeman sensitive spectral lines. On the disk, the bright photospheric and 
chromospheric emission generally overwhelms the coronal contribution, so Zeeman studies 
have been carried out above the limb, where line-of-sight integration may complicate the 
interpretation. Since Zeeman sensitivity increases at longer wavelengths, many studies in 
recent years have focussed on the promis i ng ap plication of infrared (IR) coronal emission 
lines (CELs) f)Judgelll998l: iLin et aLlDOOOl. 120041 1. Although direct comparison b etween the 
meas urement of CELs and coronal magnetic field models has been attempted flLiu fc Lin 


20081 ). the reliability of both the measurement and the model remain to be validated. 


It has been recognized that the use of microwave imaging spectroscopy can provide 


























3 


prominent diagnostics of the coronal magnetic field through the spec tral analysis of two ra¬ 
diative mechanisms : gyroresonance emissio n and free-free emission fjCarv &: Hurford 1994 : 
Ryabov et alJ Il999l: iGrebinskii et al.ll2000[) . During non-flaring times in the lower corona 


above active regions, where magnetic helds are strong (typically in the range 100 to 2000 
G), gyroresonance emission dominates in the 1-20 GHz range of microwave frequencies. The 
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ever, so far there remains no comprehensive test of the reliability of this diagnostic method 
for use on a daily basis in practical observations. With the recent advances in radio inter¬ 
ferometric instruments that will soon be capable of combining adequate spatial and spectral 
resolution for coronal magnetic held diagnostics, we provide a systematic evaluation of this 
approach, using the Expanded Owens Valley Solar Array (EOVSA) instrument prohle as 
a specihc example for quantitative comparison of the folded parameters with the original 
model. 


2. Modeling Background 
2.1. The Origin of the Mok Model 


To present a realistic challen ge in both t he im aging and spectral domains, we exploit 
a physically plausible model from iMok et al.l (1200511 . which includes a spatial data cube of 
electron density Ue, electron temperature Te, and vector magnetic held B. Galculation of 
these parameters follows two major steps: (1) the potential held model is hrst applied to 
extrapolate the magnetic held specihcation from the vector magnetogram data of one selected 
active region; (2) after the held line conhguration is determined, the thermal structure of 
the corona is computed self-consistently using a steady state heating model in which the 
volumetric heating rate is simplihed as directly propo rtional to the static magnetic held 
strength. The original data used by iMok et al.l (1200511 was based on active region 7986, 
a spatially-dispersed, magnetically weak active region observed during solar minimum in 
August 1996. To simulate a more typical active region in the corona, we reduced the spatial 
pixel dimension from 6.2 x 6.8" to 2.4 x 2.4", and scaled the total held of view (FOV) from 
Mok’s model to be 172" x 127", while keeping the original scale height. We also scale the 
magnetic held, whose original longitudinal component ranges from —221 G to 179 G, by a 
factor of 10, and reduce the electron density by 33%, to match typical ranges in the corona 





















































































above strong active regions. Althongh these changes invalidate the physical correctness of 
the active region model, for our purposes this is unimportant. We require only that the 
model contains sufficiently realistic complexity to be a fair test of the ability of extracting 
the coronal plasma parameters of the model from the simulated multifrequency radio images. 
Although the spatial resolution of the model is too low to reveal distinct loop-like features, 
it contains a physical resemblance to the actual solar corona, and for our purpose provides 
sufficient complexity for a practical evaluation. Figure [1] shows the 2D distribution of the 
longitudinal magnetic field, the electron density, and the temperature, including the cross- 
section view at the b ase of the corona, and at the horizontal center of the model. We follow 
the nomenclature in I Gary et ahl (120131) and refer to the direct images calculated from the 
model as ’’unfolded” images; and those obtained after sampling with the spatial response 
of EOVSA as ’’folded” images, ie. folded through the instrument response function. We 
continue to discuss the use of the model to generate the images at multiple microwave 
frequencies in section 3. 


2.2. Emission Mechanisms and Radiation Transfer 


We now briefly review two dominant mechanisms, gyroresonance and free-free emission, 
which have been implemented in our source code to simulate multifrequency images (a data 
cube) from the original model. Gyroresonance emission depends strongly on the plasma 
temperature (and, more broadly, the distribution function type), magnetic field strength and 
direction, and the low harmonics of gyrofrequency ub ~ 2.8x Hz are known to match the 
range of microwave frequency (1-20 GHz) in the corona where the magnetic field streng t hs can 
reach up to 2500 Gauss. Although the theory developed by iFleishman fc Kuznetsov! (120141) 
allows non-Maxwellian distributions, in this paper we consider the Maxwellian distribution 


only. For this case 


(e.g., IZheleznvakovI (llOTOh b 


he ab sorption coefficient due to a given gyro harmonic s has the form 
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( 1 ) 


from which the emissivity can be straightforwardly found using the Kirchhoff law 
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where B is the total magnetic held strength, e is the electron charge, rrie is the electron mass, 
Ufj is the refractive index of the wave-mode a, 9 is the angle of the held to the line of sight. 
To- and To are the compone nts of the wave polarizati o n vector (al l variable dehnitions follow 
the conventions adopted in iFleishman fc Kuznetsov! (120101.12014J) ). Note that the presence 
of the exponential factor in Eq ([1]) makes the gyro opacity exponentially small everywhere 
beyond a very narrow range of frequencies very close to the gyro harmonics svb- For emission 
produced at a given frequency in the solar atmosphere, where the magnetic held varies with 
height, this narrowness translates to a narrow range of heights, called a gyrolayer, which 
makes a dominant contribution to the opacity at this frequency. For typical conditions of 
the corona, where the magnetic held obeys the resonant condition 


eB 


V = SVb = S 


2'Kmf.c 


2.8 X IOST Hz, 


(3) 


the gyrolayers are only about 100 km thick. For the solar corona, s is a small integer, typicall y 


1, 2 or 3, although in some cases s = 4 may have appreciable opacity fjWhite fc Kundu 


1997|). 


Wheneve r the gyroresonance emission is small, the free-free emission makes a dominant 
contribution flGarv fc Kellerl 12004 ). In a magnetized plasma the free-free emission of a given 
eigen-mode has the form (Zlotnik 1968a l^: 
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where nu is the total number density of all singly ionized atoms (mainly-protons)l§, nn 
He — and UBeiii IS the number density of the fully ionized Helium 
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V = sin"^ 9 + 4m(1 — vY cos^ 9^ 
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( 6 ) 
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^Expression of the factor respo nsible (along with in the d enominator) for polarization of the free- 
free emission is given in textbook by Fleishman fc Tootygin ( 201 ,*ll) but with a typo resulted in the wrong 
general sign of the expression. 


^If no independent information about ionization states (e.g., non-LTE) is supplied by the given model, 
the code computes the ionization states of hydrogen and helium using the Saha equation based on the 
electron density and temperature input. Note that for the coronal temperatures and coronal abundances, 
Ue(n// +^nHeiii) ~ l.llng. 


































6 


z/pe = is the electron plasma frequency. For X-mode, a = —1; for 0-mode, 

a = +1, Ac is the Coulomb logarithm. 


Unlike the case of gyroresonance opacity, the free-free opacity is generated over a much 
larger LOS, although when seen optically thick against the solar disk it is heavily weighted 
toward the lowest, densest part of the ray path; the lo wer the frequency, the higher the 


effective formation height (e.g., iLonkitcheva et al.l 120151 ) of the free-free emission ranging 


from chromosphere at high frequencies (> 10 GHz) to corona at lower frequencies. 

In principle, the resulting gyroresonance and free-free optical depth and emerging radia¬ 
tion can be calculated for each volume element (voxel) from any model that provides electron 
number density, electron temperature, and vector magnetic held in each voxel. However, the 
size of the voxel height is typically much larger than the dept h of a gyro layer. For that rea- 
son, the computation engine (see the complete description in iFleishman fc Kuznetsov! 120141) 
linearly interpolates the magnetic held vector between the neighboring voxels to precisely 
locate the positions of all gyro layers of interest (the default maximum value Smax = 7 is 
adopted but can easily be adjusted by the user) to calculate the absorption coefficient (0]), 
and then solves the equation of radiation trans fer, taking account of frequency-dependent 
mode coupling as described in iNita et al.l (120141) . The equation of radiation transfer solves 
for the intensities of the X and O modes, which have elliptical (i.e., not necessarily circular) 
polarization when they have just left the data cube. We take into account, however, that 
the waves continue to propagate through the corona with declining electron density and 
magnetic held and so the elliptical polarization necessarily evolves t owards a truly circular 
one d ue to the ehect of ’limiting polarization’ fjZheleznvakovi Il970l: iFleishman fc Toptvgin 
20131 ): thus, we identify the intensities of the X and O modes with the circularly polarized 
waves observed at the Earth. 


3. Modeling Scheme 
3.1. Generation of Unfolded Images 

In order to simulate the multifrequency unfolded images, we calculate the radio emis¬ 
sion for each LOS pixel, at 64 logarithmically-spaced frequencies in the range 1-18 GHz. By 
considering each voxel as a homogeneous source, the model parameters at the center of a 
given voxel can be used to calculate the resulting, frequency-dependent total optical depth 
A,a = 'Tgr,crATff,a ffi two polarizatious at that voxel. Subsequently, the total brightness tem¬ 
perature of a given pixel at each frequency can be obtained through appropriate integration 
of radiative transfer in the corresponding voxels along the LOS. These calculations are rou- 
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tinely done within the GX_siniulator tool flNita et alJl2014i) to produce a full set of simulated 
2D brightness temperature maps at the chosen frequencies. The hrst column from left of 
Figure [2] shows a sample of 6 of these RCP (unfolded) images out of the 64-frequency data 
cube. The similar LCP (unfolded) images are also shown in the third column of Figure |2l 


3.2. Generation of the Folded Images 


In an actual observation of gyroresonance emission, from which the coronal magnetic 
held is to be inferred, the target active region brightness is competing with that of other 
active regions and the radio-bright solar disk, so that the uv visibilities contain contributions 
from this complex spatial structure. This can contribute to confusion in the reconstructed 
image due to hnite sampling in the uv plane. This problem is especially true in EOVSA’s 
case, since it consists of a limited number of small antennas (a total of 13 2.1-m antennas). 
At 18 GHz, the primary beam size of one EOVSA antenna is about 33 arcmin, increasing to 
18 times larger at 1 GHz following the diffraction limit. Therefore, to provide a fair test of 
the instrument performance, we embed the unfolded images, without loss of generality, into 
the central region of a full-disk model. Figure [3] is one of the 128 simulated full disk images 
(64 frequencies and 2 polarizations) at one selected frequency and polarization. 

The frequency-dependent full-disk model is generated as follows. As a base map we use 
SOHO/EIT EUV images at 171 and 195 A. The ratio of these two images is converted to 
emission measure and temperature (in the range 0.6-1.7 MK) using standard EIT software 
and assuming a coronal abundance for Fe. The date of the EUV images is 1999 Apr 11 
(at 19 UT), chosen because it was available from an earlier (unpublished) study. From the 
emission-measure and temperature maps, we can generate a model radio image (in units 
of brightness temperature for convenience) from bremsstrahlung that correctly accounts for 
optical depth variation with frequency (i.e., at low frequencies active regions are optically 
thick at the temperature of the model). 


The EIT model only reproduces bremsstrahlung from coronal plasma. The radio im¬ 
ages also see a background from the cooler lower atmosphere that we add as a frequency- 
dependent disk of a brightness temperature that results from hts to the solar radio spectrum 
at solar minimum, assuming an effective disk radius determined f rom hts to t he size of full- 
disk images at frequencies where such images are available (e.g. lGarvlll996L — at 4.5 GHz 
the radio limb is 30" above the photosphere). 


Gyroresonance sources typically he above strong magnetic helds in the photosphere. In 
order to simulate them we use an MDI line-of-sight magnetogram, and apply a frequency- 








dependent threshold in held strength as well as a dilation transformation that gives gyrores- 
onance sources of about the same size as are seen in the actual radio images on 1999 April 
11. The gyroresonance sources are assumed to be optically thick at 2.5 MK. As an ansatz 
to generate plausible polarization models, gyroresonance sources with positive polarity are 
added to the right circular polarization model (replacing the other sources, since the gyrores¬ 
onance sources are optically thick), while those with negative polarity are added to the left 
circular polarization model, i.e., all gyroresonance sources are assumed to be 100% circularly 
polarized, which is not true of actual sources. Lastly, for each disk brightness model with 
suitably realistic spatial complexity, the unfolded images from the active region model at the 
corresponding frequency are added at a spatial location near disk center, as shown for 4.15 
GHz in Figure [3l 


To produce the instrument-folded image datacube, we use EOVSA as a test case, i.e., 
a 13-antenna array (78 baselines) with imaging from 1 to 18 GHz. Figure 0] shows the uv 
coverage of EOVSA at 18 GHz, used for a simulation covering 12 hours of observation on a 
date whe n the Sun is at -|- 15° declination. Then, an automatic procedure using the Miriad 
package fjSault et al.lll995l) is followed: (1) the brightness image of the disk with embedded 
active region model is Fourier-transformed into the uv plane to generate the resulting model 
visibility function V(M,n); (2) the true visibility function is sampled with the 78 baselines 
of instrumental prohle (the samples at each frequency and time V(u, v, t, v) are termed 
visibilities), shown in FigurelH Since the effect of bandwidth smearing can potentially reduce 
the quality of reconstructed images, we include a bandwidth of 50 MHz in the uv sampling, 
which slightly exceeds the maximum limit of EOVSA’s bandwidth (variable up to 40 MHz); 
(3) thermal noise from the system, 4500 K, is added to the visibilities to simulate the effect 
of the bright solar emission on system temperature, and (4) the image is reconstructed using 
Miriad's standard Maximum Entropy Method (MEM) routine maxen. By repeating (l)-(4) 
for each frequency and polarization, a full set of RGP and LGP images is generated at all 
sampled frequencies. 


The second and fourth columns of Figure [2], shown at the same frequencies as in the 
unfolded images in the hrst and third columns correspond to the relevant portions of the 
simulated (folded) images in RGP and LGP, respectively, after being sampled with the uv- 
coverage of the EOVSA instrument (Figure 0]) . In general, the features of the active region 
from the folded maps follow the unfolded ones very well, especially in higher frequency 
range, say above 4 GHz. However, the folded images tend to be lacking in details at lower 
frequencies: around 1-2 GHz, the original sharp drops in brightness temperature within the 
opacity holes directly over the sunspot areas are fairly smooth in the folded maps in both 
polarizations. This is due to the poorer spatial resolution at lower frequencies, as set by the 
maximal projected baseline i?A,proj of the antenna arrays. We show in the following that 
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even with limited visibility coverage, EOVSA can provide reliable spectra for the coronal 
magnetic held measurement in such a complex active region structure. 


4. Results and Analysis 


4.1. Comparison of the Brightness Temperature Map and Spectra 

We start with the elementary radiative transfer problem of radio emission and extend 
our discussion over emission features from our active region model in the unfolded Ti, maps. 
In microwaves, the intensity of the source can be approximated by the Rayleigh-Jeans law 
in terms of ly ~ 2ku‘^Tb/c^, where Tb is the brightness temperature, which is equivalent to 
the temperature of a blackbody having the same brightness. The radiative transfer equation 
in terms of brightness temperature is given by: 



( 8 ) 


where r is the maximum opacity along a given LOS. The second term represents the con¬ 
tribution of a background such as the optically-thick chromosphere. For an optical-thick 
thermal source under thermal equilibrium (r 3> 1), the observed brightness temperature is 
approximately equivalent to the local coronal electron temperature Tb ^ Tg. For gyroreso- 
nance emission alone, which is the dominant mechanism in the strong magnetic held (above 
100 Gauss), emission is conhned to the harmonic layers satisfying equation (1). Therefore 
the measured brightness temperature along a given LOS is contributed by the outermost 
gyroresonance layer that is optically thick at a given observing frequency. 

It is instructive to view the unfolded data cube vertically to examine the microwave 
spectral properties at particular LOS positions, from which the physical parameters can in 
principle be extracted through spectral analysis. Figure [5](a)-(b) show the unfolded and 
folded images in RCP at 4.3 GHz, which corresponds to the magnetic held strength |i?| ~ 
512 G of the third harmonic layer. The numbered points indicate the locations of six sampled 
spectra in Figure [5](c)-(h). Red and blue represent the right- and left-circular polarized (RGP 
and LGP) spectra. The unfolded spectra (solid lines) in the two polarizations at a given LOS 
diher due to diherences in the highest optically-thick harmonic, which can be read directly 
from the unfolded spectra. 

To aid in further discussion, we show in Figure [H] the parameters of the model as a 
function of height along the 6 lines of sight marked in Figure [5l The inferred brightness 
temperature of point 1 (red curve), sampled near the loop-top region, is around Tr ^ 2.3 x 10® 
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K at 4.3 GHz (Figure (St). The height-dependent magnetic held of the model (Figure [6 )d) 
shows that the corresponding 3rd-harmonic held strength \B\ = 512 G (thin dotted lines in 
Figure [Hb) meets the red, point 1 curve at a coronal height of ~ 38 Mm, where the coronal 
electron temperature is Tg ~ 2.7 MK (Figure [ntL)- At point 3, both RGP and LGP spectra 
(Figure [5^) show steep drops in brightness temperature at 6 GHz. The LGP spectrum then 
maintains a moderate brightness above 6 GHz, with a second drop in brightness around 9 
GHz. This corresponds to the next higher (s = 3) harmonic layer, which is partially optically 
thin, dropping out of the corona. The 2:3 ratio of the frequencies of the two sharp drops 
in the spectrum indicates that the turn-over frequency of 6 GHz corresponds to the second 
harmonic layer. Accordingly, at 6 GHz the held strength derived from s = 2 is |H| ~ 1071 G, 
which is the limiting held strength along that LOS (Figure [6 ]d), and occurs at a model height 
of 77 ~ 1 Mm and temperature Te ~ 10^ K. This accounts for the sudden drop in brightness 
at that frequency. Away from the active region (points 4-6) where the 512 G isogauss layer 
does not he in the corona (the curves for points 4-6 all fall below the thin dotted line in 
Figure [6]d), free-free emission is the dominant mechanism. For example, at point 4 above 
3 GHz (Figure O), the gradually-decaying spectrum suggests that bremsstrahlung is the 
primary emission mechanism. 


It is also important to interpret the nature of the sudden depression in the brightness 
temperature over the sunspot regions in both polarizations (e.g Tr at point 2 and 3, Fig¬ 
ure [5^) in terms of local spectral behavior (Figs. [5li,e). The reason for the depression is 
that when the angle 9 between the vector magnetic held and the LOS approaches 0° or 
180°, the corresponding gyroresonance layer becomes optically thin. This can occur in the 
0-mode emission at larger angle 9 = 9c than for X-mode emission (jWhite fc Kundul 119971: 


Gary fc Kelleiil2004j ). At point 2, the angle 9 along the LOS (magenta line in Figure [6t) is 


less than 20°, which falls below the critical angle 9c in the O mode (Tr). Thus, there is no 
signihcant brightness contribution to Tr from the optically thin harmonic layers along this 
LOS, leaving only free-free emission in the Tr spectra (Figure Oi). In comparison, point 3 
(Figure [5^) also shows similar spectral behavior in the 1-2.5 GHz range, except for a high 
brightness temperature bump that peaks at 5.2 GHz. We showed earlier that this corre¬ 
sponds to second-harmonic emission, which yields a held strength B = 928 G (thick dotted 
line in Figure [nt*)- The corresponding angle in Figure [Ht shows 9 = 33.2°, which is at least 
10° larger than at point 2. In this range, 9 approaches 9c, increasing the opacity and hence 
the brightness, as frequency increases toward 6 GHz. The variation of brightness tempera¬ 
ture with viewing angle is rich in diagnostic power, but does complicate the interpretation 
of spectra (see section 4.3). 


Finally, the dashed lines in Figure [5t-f show the effect on the spectrum of folding the 
model images through the instrument prohle. Some folded spectra match the unfolded spec- 
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tra (solid lines), but due to the limited uv coverage and the resulting hnite spatial resolution, 
the folded spectra tend to be both more extended and lacking in spectral details. Excellent 
agreement can be found in the center of the active region, where the spatial brightness tem¬ 
perature gradient is relatively smooth. For example, the folded spectra from point 1 show 
excellent agreement in both polarizations over a broad frequency range. In contrast, the Tr 
or Tl spectra from point 2 to point 6 in Figure [5li-h show varying degrees of distortion. At 
turn-over frequencies, the dominant polarization (that due to the higher harmonic) of the 
folded spectra generally agree with the unfolded spectra down to 25% ( 0.5 MK) of the peak 
brightness ( 2 MK) for all except point 5 and 6, while the other polarization is good only 
up to about 50% of the peak. This suggests that (not surprisingly) robust estimation of the 
coronal magnetic held strength may be limited to the brighter parts of the radio images, 
down to some limiting brightness temperature. We discuss this further in the next section. 
Because the spatial resolution continues to decrease below 3 GHz, spectra with sharp varia¬ 
tions in frequency there tend to be distorted (e.g. the spectra of point 2 and 3 in RCP, and 
point 6 in both polarizations). For instruments with more antennas and longer baselines, 
such as the Jansky Very Large Array (VLA), the Chinese Solar Radioheliograph (CSRH), 
the upgraded Siberian Solar Radio Telescope (SSRT), or the future Frequency Agile Solar 
Radiotelescope (FASR), the observed spectra will more-closely approach the quality of the 
unfolded spectra of Figure O and hence a spectrally-based analysis will be possible. We 
discuss this further in section 14.31 


4.2. Coronal Magnetography using an Image-Based (Level-0) Method 


As described in section ITTl the gyroresonance spectrum is characterized by a sharp drop 
in brightness temperature from coronal to chromospheric values. The spectral shape is merely 
a reflection of the fact that the temperature structure of the solar atmosphere along any LOS 
has a steep drop that occurs at the transition region, where the frequency-dependent gy¬ 
roresonance layer drops out of the corona. Therefore, the magnetic held strength close to the 
base of the solar corona can be determined from the relevant harmonic of the gyrofrequency 
through equation (1), for lines-of-sight where gyroresonance emission dominates. Such a 
technique has been attempted using observational data to provide a rough estimate of mag¬ 
netic helds at a few discrete levels co rresponding to the hxed frequencies avai lable for imaging 
in previous studies flLee et al.lll998l). and at a larger number of frequencies flGarv fc Hurford 
1994JdTun et al .11201 It iKorzhavin et al.ll201Cll . However, as Figure E] shows, for a realizable in¬ 
strument with hnite spatial resolution the observations compare best with the model-derived 
images only for the higher brightness temperatures (say above about 25% of the peak bright¬ 
ness in the case of EOVSA). This suggests use of an image-based approach. For simplicity. 
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we refer to image-based analysis as the Level-0 method. In this section, we hrst introdnce 
the principle of the Level-0 method, and apply it to the simnlated EOVSA data. Then, we 
systematically evalnate its performance in dedncing the magnetic held strength. 

The Level-0 method uses three steps to deduce the coronal magnetogram, illustrated for 
three frequencies in Figure [71 (i) for each radio map at a given frequency, the contour of a 
chosen brightness temperature is used to represent the “edge” of the gyroresonance surface 
of the presumed harmonic value s = 3; (ii) the resulting magnetic held strength along this 
temperature contour can be deduced from the observing frequency used to produce the image, 
via equation (1), assuming s = 3; (iii) a set of magnetic held contours, derived in the same 
way from multiple-frequency radio maps, can be overlaid to construct a hnal magnetogram at 
that chosen temperature level. Note that for each LOS, the brightness temperatures in RCP 
and LCP usually diher due to the dominant contributions being from diherent harmonics. 
To resolve this temperature ambiguity, we choose the dominant polarization, i.e. the one 
with the higher-frequency cut-oh at the chosen threshold temperature, to correspond to the 
X mode. Generally, the choice of harmonic s = 3 is expected to be correct for the X mode 
in coronal conditions in active regions, except over sunspots when the angle of the magnetic 
held to the LOS is close to 0° or 180° (this is usually only a signihcant factor for sunspots 
near disk center, as in our example). In these regions, the lower optical depth causes the 
harmonic ratio to be 2:1. Thus, we begin by making the usually correct assumption that 
s = 3 for the X mode, keeping in mind that in regions near the cores of sunspots away 
from the limb the correct choice may be s = 2. Note that our assumption of s = 3 as 
the highest optically-thick harmonic is the correct choice at points 1, and 4-6 in Figure O 
but is not correct at points 2 and 3 where the harmonic ratio is 2:1 based on the unfolded 
spectra. Nevertheless, following the above steps. Figs. [7^-c show the folded images in X 
mode polarization at three frequencies. A relatively high threshold temperature of 0.72 MK 
has been chosen to characterize the edge of the gyroresonance layer, as indicated by the 
dashed contour over each hgure. Figure [Tli shows the corresponding magnetic held strengths 
along these overlaid contours. By applying the same technique for the many images of the 
data cube over the entire frequency range, we obtain the overall magnetogram shown in 
Figure 03. 

We can evaluate the reliability of the Level-0 approach by comparing the Level-O-derived 
coronal held map, Sfoid in Figure |8 ]d, with the actual magnetic held strength in the model, 
Bmodei, at the height where the temperature Tg = 0.72 MK occurs, shown in Figure [8^. The 
Bfoid map extends over only the central part of the region, because the relatively high thresh¬ 
old temperature does not extend to the periphery. A dashed contour at 600 G, obtained 
from Rmodei, is plotted over each hgure for comparison. To provide a more quantitative com¬ 
parison for each LOS, we calculate the percent error (i?foid —7^modei)/7^modei- Figure [HI:; shows 
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the resulting error distribution. Each color represents a certain error range, as indicated by 
the color bar. Above ~ 120 G, both ilmodei and i?foid qualitatively agree very well over the 
active region where an error tolerance of < 30% can be widely achieved. Note that abrupt 
variations from positive to large negative error, indicated by the blue and purple segments, 
occur in the regions near the sunspots, where the magnetic fields are nearly vertical (see 
point 2 in Figure [61::), so that the harmonic ratio becomes 2:1 and our assumption of s = 3 is 
violated. For quantitative comparison. Figure [Hll shows the Sfoid (symbols) and Smodei (solid 
line) sampled in a horizontal cut across the central part of the active region, as indicated 
by the dashed line in Figure [HI;. The prohle of Ufoid, represented by asterisks, hts very well 
with the i^modeh except for the clearly discrepant, underestimated magnetic fields around 
the sunspots. A subjective correction of the discrepant points by a factor of 3/2, which 
comes from using a lower harmonic (s = 2) in equation (1), shows improved agreement as 
indicated by open circles. In general, the result from the simple Level-0 procedure can fully 
exploit the un-coverage of EOVSA to generate a reasonably reliable coronal magnetogram 
despite the limited spatial resolution of the folded data, although it is to be expected that 
there will be localized deviations near sunspots where the presumption of harmonic s = 3 is 
inappropriate. 

Although in principle we can derive a coronal magnetogram at any temperature level, 
we argued earlier based on Figure [5] that it is most favorable to choose a higher probing 
temperature to reduce the spectral distortion caused by two factors: the spatial resolution 
and the free-free emission, both of which can induce a gradual tail in the folded spectra above 
the turn-over frequency. It follows from equation (3) that the free-free emission depends on 
the square of the electron density Ug. As an illustration, we reran the model after increasing 
He everywhere by a factor of 3 (i.e., an order of magnitude increase in free-free emission), 
with results shown in Figs. [8^-f. In this case, the magnetic field Sfoid derived from the Level 
0 method deviates significantly below ~ 600 G, corresponding to a turn-over frequency of 5 
GHz at the third harmonic. This is due to much of the lower-held region being covered up by 
overlying free-free emission, exacerbated by the lower spatial resolution at lower frequencies. 
Finally, since the Level 0 magnetogram is derived at a constant temperature level, the 
downside of using a higher value for this threshold is that it corresponds to greater and more 
variable heights in the corona. Observations with better spatial resolution (e.g. VLA, GSRH 
or FASR) will support the choice of lower threshold temperatures (although in that case the 
spectrally-based “Level-1” approach may be better, as discussed in the next section). 

Going back now to our original model with unenhanced density. Table 1 lists key quan¬ 
tities at our illustrative locations 1-6 in Figure [Hl(c). Golumn (6) shows the percent error 
of magnetic held derived from the folded images, Rfoid, from Level 0 listed in column (4) 
with respect to the model, Rmodei, in column (2). The total magnetic held at point 1 shows 
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remarkable agreement because the turn-over frequency i/foid agrees well with z/unfow at the 
brightness temperature Q.72 MK. At points 5 and 6, there are severe losses of spatial and 
spectral resolution in the low frequency range of the folded spectra, which results in a higher 
turn-over frequency and overestimated magnetic held strength -Bfoid- The B values at points 
2 and 3 deviate from the model by —34%, and —16%, respectively, due to the use of an incor¬ 
rect harmonic number. In real observations, an empirical correction to the choice of harmonic 
might be possible, but any automated approach will require experience with the actual data. 
In the next section, we further discuss the potential for more advanced, spectrally-based 
analysis, which we refer to as the Level-1 method, to correct for misdiagnosed harmonics. 


4.3. Coronal Magnetography using a Spectrum-Based (Level-1) Method 

We now briehy explore the prospects for obtaining more accurate coronal magnetograms 
through Level-1 analysis, which makes use of the spectra along each LOS in the two indepen¬ 
dent circular polarizations. Although EOVSA will not have the spatial resolution required 
to use this method in most cases, the VLA and some future instruments will. Along each 
LOS, the two distinct propagating electromagnetic modes, the X mode and O mode, can 
be utilized to remove ambiguity in the appropriate harmonic. As noted in section 4.1, it is 
commonly seen that the local spectra in O mode fall off at a lower frequency than in the X 
mode due to the dominance of different harmonics for the two modes. The Level 1 method 
can automatically detect a correct harmonic for each LOS in those cases where the steepest 
drops of the two polarized spectra occur, and use the corresponding turn-over frequencies 
vr and Vl to match particular harmonic ratios (2:1 or 3:2, and possibly 4:3 in some cases). 
In general circumstances where vr k. ul, a harmonic ratio of 3:3 is likely since the third 
harmonic then may be optically thick in both polarizations (e.g., when the viewing angle 
is nearly perpendicular to B). We test this approach to resolving the harmonic number in 
the absence of frequency-dependent resolution by applying our algorithm to the unfolded 
data cube. Figs. |9t-b show the two examples of local spectra (obtained at LOS points 3 
and 5) where the correct harmonic ratios, 3:2 and 4:3 respectively, have been successfully 
determined. Figure IH:: shows a percent error map, in which the Level 1 coronal magne¬ 
togram deduced from the unfolded image cube is compared with the Bmodei map. Columns 
(3) and (7) in Table 1 represent the unfolded magnetic field derived from the Level-1 method, 
and the resulting percent error. Figure [Hll shows the corresponding harmonic map detected 
by the algorithm. It has been discussed in section 4.1 that O mode might become totally 
optically thin at all observable frequencies (e.g. point 2), rendering the harmonic determi¬ 
nation via the O mode emission unusable. In Figure [Oji, this corresponds to the cores of 
the two sunspots, which are still presumed to have s = 3 as the appropriate harmonic but 
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are surrounded by the second-harmonic pixels. Therefore, it is reasonable to replace them 
with s = 2. Point 2 has been manually corrected using s = 2, yielding a percent error 
that reduces to —3% (column (3) in Table 1). In general, the Level-1 method can resolve 
the misdiagnosed harmonics over the entire map, especially within the central active region. 
Figure [H:: shows that with higher resolution, systematic effects due to finite resolution can 
be significantly reduced to provide a more accurate result over a much larger area in the 
active region compared to Figure [HI::. Away from the active region, where B is lower than 
a few hundred gauss, the unfolded magnetogram can still provide 30% error on average. 
Since this algorithm relies heavily on the steepness of the polarized spectrum, it occasionally 
fails to diagnose higher harmonic ratios. For example, at point 5, the steepest slopes of 
both polarizations are detected at 2.3 GHz and 2.9 GHz (Figure [9 ]d), whose ratio is close to 
s — 1 : s = 4 : 5. Instead, a harmonic ratio of 3:4 is clearly the correct one because the steep 
drop of Tr at 1.5 GHz should be interpreted as the second harmonic. In brief, a more pow¬ 
erful algorithm could be developed to improve the Level 1 result as higher quality radio data 
become available in the future. W e also anticipate th at a forward fitting approach similar to 
that described for flaring loops in [Gary et ah! (120131) can be developed, which will have the 
advantage of providing additional coronal parameters from the gyroresonance spectra, such 
as density and angle of B from the LOS. Such spectral htting can also be used to explore 
whether alternative, no n-Maxwellian particle distribu tions such as k— and n— distributions 
exist in active regions ([Fleishman &: Kuznetsovll2014f) . 


5. Discussion and Conclusions 

In this paper, we have demonstrated that EOVSA will be able to measure coronal 
magnetic field strengths in active regions, in areas where gyroresonance emission dominates. 
To evaluate the reliability of recovering the held strength, we performed the following steps: 
(1) We adopted an active region model, in which 3D coronal temperature, electron density 
and vector magnetic held are complex enough to provide an adequate challenge to test the 
feasibility of deriving magnetic held strength from microwave interferometric imaging. (2) 
We simulated a full set of unfolded images at multiple radio frequencies from 1-18 GHz, 
in which two dominating mechanisms, free-free emission and gyroresonance emission, have 
been included. (3) To provide a fair test of the resolving power of EOVSA, we embedded 
the active region into a realistic, frequency-dependent solar disk model. (4) A standard 
procedure was followed to generate two sets of model visibilities—those directly from the 
images (“unfolded”) and those folded through EOVSA’s instrumental prohle after 12 hours 
of rotation synthesis (“folded”). Sets of folded images were then reconstructed from the 
maximum entropy algorithm implemented in Miriad. 








We have used these reconstructed data cubes to investigate two alternative approaches 
to deriving coronal magnetic held maps, the simple Level-0 (image-based) method, and 
the more sophisticated Level-1 (spectrum-based) method. We have shown that even with 
limited uv coverage, in our case 78 baselines over a 12-hour observation, the Level-0 method 
can provide a reliable coronal magnetogram covering most areas within the active region. 
However, the deduced magnetic held strengths showed deviations from the model in a few 
areas near the cores of sunspots where the true harmonic in the X mode was not the presumed 
s = 3 value. To resolve the harmonic ambiguity, higher resolution radio data are required, 
although some ad hoc correction is possible simply by making the s = 2 assumption in these 
sunspot areas (although this is inhuenced by the fact that we put the model active region 
near disk center so that the magnetic held over the umbrae is nearly parallel to the line of 
sight). We have briehy explored the power of the Level-1 method using direct model images 
of the unfolded data cube. The principle of this technique relies on a simple interpretation 
of the faithfully reconstructed polarized spectra, in which the s : s — 1 harmonic frequency 
ratio can be determined in order to identify the relevant harmonic layer. The result has 
shown that most misdiagnosed harmonics have been corrected with signihcantly improved 
accuracy in a wider area. 

The spectrum-based Level-1 method requires excellent spatial resolution as well as high 
image hdelity. With the advent of new radio facilities capable of multi-frequency solar radio 
imaging, such as EOVSA and the Jansky VLA, as well as future solar-dedicated instruments 
such as CSRH, the upgraded SSRT and FASR, the solar community is on the verge of 
having true coronal magnetograms that will grow increasingly reliable as analysis techniques 
are rehned through experience. We expect that breakthroughs on long-standing solar topics, 
like the coronal structure of active regions, coronal heating and magnetic reconnection, can 
be achieved in the near future. 

A hnal word about the interpretation of such coronal magnetograms is in order. As for 
any remote sensing method for measuring magnetic helds, the magnetic held is measured 
at the location of emission formation rather than at a given height. For photospheric lines, 
the range of heights is small and is generally ignored, but height variations become more 
important for chromospheric and coronal techniques such as the one described here. As we 
have seen, the spatially-resolved gyroresonance radio spectrum provides the measurement 
of Ti,{u), which by virtue of the resonance [u ^ ub ^ B) and optically-thick conditions 
(Tfe —)■ T) becomes T{B) (or equivalently B{T)). By making use of the sharp drop in 
temperature that occurs at the base of the corona, we can constrain the height somewhat, 
but the precise height to which the measurement applies can vary by a few thousand km. 
Therefore, coronal magnetograms will never be as straightforward to interpret as those for 
the photosphere, and their use will be most helpful in combination with 3D models. This 
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paper has shown that B{T) can be reliably measured using radio interferometry data of solar 
active regions, which is the essential starting point for fruitful comparisons with 3D models. 
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1250374, and NASA grants NNX11AB49G and NNX14AG87G to New Jersey Institute of 
Technology. 
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Table 1: Comparison of derived parameters for the model, unfolded, and folded spectra, 
using Level 0 and Level 1. See detail description of each parameter in section 4.2-4.3 


(1) 

Point 

(2) 

Bmodel (G) 

(3) 

Bun fold (Gl) 

(4) 

Bfoid (G) 

(5) 

s 

(6) 

AS", 

(7) 

AB4 

1 

1476 

1422 

1387 

3 

-6.0 

-3.7 

2 

1675 

1667 

1132 

(2) 

-33.7 

-2.6 

3 

997 

1141 

838 

2 

-16.0 

14.4 

4 

394 

411 

385 

3 

-2.1 

4.3 

5 

208 

265 

223 

3 

7.4 

27.5 

6 

134 

152 

160 

3 

19.1 

12.8 


Note. — The actual harmonics in column (5) are obtained from the Level 1 algorithm, except point 2 
whose harmonic is manually corrected to minimize the percent error. 
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Fig. 1.— (a)-(c) The longitudinal magnetic field electron density N^. and temperature 
Tf. distribution at the model base corresponding to the transition region. The field of view is 
172" X 127". (d)-(f) The corresponding height-dependent distribution of coronal parameters 
along the horizontal line in (a)-(c). 
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Fig. 2.— A comparison of the unfolded and folded images in both polarizations at six 
selected frequencies from 1-18 GHz. Images from left to right: the hrst and the third column 
are the unfolded images, which are derived from the radio emissions of the Mok’s model, 
in RCP and LCP respectively; the second and the fourth columns are the corresponding 
restored images, after folding through EOVSA. The total FOV of each map is 172" in the x 
direction, and 127" in the y direction, with uniform spatial sampling of 2.4". 
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Fig. 3.— A full-disk unfolded radio image in RCP, calculated for 4.15 GHz. The active 
region model is embedded near the center of the disk. 
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Fig. 4.— The u,v coverage (EOVSA sampling function) at 18 GHz, which uses an hour angle 
range -6 hr <h <6 hr and +15° declination to produce the folded images. 
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Fig. 5.— (a)-(b) The RCP image at 4.34 GHz, before and after folding throngh the EOVSA 
instrument, respectively, (c)-(h) Comparison of the sampled spectra at six locations. The 
spectral curves in red and blue represent RCP and LCP, respectively. The solid lines, and 
dash lines with symbols represent the unfolded and folded spectra respectively. 
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Fig. 6.— Variations along the 6 lines of sight throngh the model of (a) the electron tem- 
peratnre Te, (b) the magnitnde of the total magnetic held \B\, and (c) its absolnte angle 
l^l as a fnnction of coronal height. The sampling points 1-6 are indicated in red, magenta, 
orange, green, cyan, and bine, respectively. The part of the height range corresponding to 
emission from 1 GHz to the tnrn-over freqnency in the Tx spectrnm of each line of sight is 
accentnated by thicker lines. The pairs of dotted lines are discnssed in the text. 
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Fig. 7.— (a)-(c) The brightness temperature maps at representative frequencies 1.20 GHz, 
4.98 GHz and 12.47 GHz, respectively. The brightness temperature at each pixel shows the 
local value of either RGP or LGP brightness temperature, whichever is higher. A contour 
of a given brightness temperature Tb ~ 0.72 MK, as represented by the dashed line, is 
plotted to outline the sharp edge of gyroresonance surface. The total magnetic helds along 
the contour are derived from equation (1). (d) The magnetic held contour map based on the 
three frequencies in 6(a)-(c) at Tb ~ 0.72 MK. Similar contours at other frequencies hll out 
the hnal Level 0 coronal magnetic held map. 
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Fig. 8.— (a) Result from the model: the magnetic held strength Bmodei evaluated at coronal 
temperature 0.72 MK. (b) Result from Level 0 method: the magnetic held map Bfoid 
at brightness temperature 0.72 MK, which is derived from the 64-frequency folded data 
cube, (c) Schematic diagram of percent error in quantitative comparison of (b) with (a), 
(d) Comparison of Bmodei (solid line) with Bfold from the Level 0 method (asterisk symbols) 
along the horizontal line in Figure 8(c). The open circles represent manual corrections using 
a harmonic s = 2. Figure 8(e) and (f) are obtained via similar procedures as Figure 8(c) and 
(d), but the electron density has been scaled to 3 time larger. The resulting strong free-free 
component reduces the accuracy widely for helds below 600 G. 
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Fig. 9.— (a)-(b) Two examples of polarized spectra from the unfolded data cube, where the 
gyrofrequencies, corresponding to the sharp drops in RCP (red) and LCP (blue), are found to 
match correct harmonic ratios, (c) The percent error map in which Bunfoidi deduced from the 
Level 1 method, and Bmodei are compared. Most misdiagnosed harmonics are automatically 
corrected within the central region, (d) The corresponding harmonic map diagnosed from 
the unfolded image datacube using the Level 1 algorithm. The black region, labeled s = — 1 
on the colorbar, shows the region where no harmonic could be undetermined, because Tf, is 
less than the probing temperature. 
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